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■ Abstract. In order to describe large transverse momentum {pj) distributions observed in high 
energy nucleus-nucleus collisions, a stochastic model in the three dimensional rapidity space is 

, introduced. The fundamental solution of the radial symmetric diffusion equation is Gaussian- 

■ like in radial rapidity. We can also derive a pr or radial rapidity distribution function, where a 
\ distribution of emission center is taken into account. It is applied to the analysis of observed large 

^-H , pj distributions of charged particles. It is shown that our model approaches to a power function of 

Pt in the high transverse momentum limit. 
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O : INTRODUCTION 
^ , 

vn ■ 

O ' In high energy nucleus-nucleus collisions, colliding energy ^JsNN per nucleon grows up 

^ . to y/Sf^M = 200 GeV, and secondary particles with transverse momentum pj more than 

Ph! 10 GeV/c are observed. The observed pr distributions have long tail compared with 

p exponential distribution in 

1^ ■ As is well known, the fundamental solution of stochastic process is Gaussian if 

'. variables in the Euclidian space are used. Therefore, as long as we consider stochastic 

I equations in the transverse momentum space, it is very hard to describe observed pr 

^ ■ distributions. This fact suggests that a relativistic approach to the stochastic process 

^ ! would be needed. 

" " ' In Ref. lUl], an empirical formula for large pr distributions at polar angle 9 = n/2, 

\B=n/2 = Aexp[-yj-/{2LT)], 



d^p 



1 E + \pj\ 

yT = o^n-- . r, (1) 

2 £ — Iprl 



was proposed from the analogy of Landau's hydrodynamical model. Polar angle 6 
(0 < < tt) is measured from the beam direction of colliding nuclei or incident particles. 
In Eq. (HJ), E denotes energy of an observed particle, Lt is a parameter, and yx is called 
the "transverse rapidity". Equation ([T]) well describes large pj distributions for p + p ^ 
-\-X and p + p 7t^ -\-X. However, it cannot be derived from the hydrodynamical 
model. 

As for the relativistic approach to the stochastic equation, we consider the diffusion 
equation in the three dimensional rapidity space or Lobachevsky space, which is non- 



Euclidean. In order to specify longitudinal and transverse expansions, it would be ap- 
propriate to consider the diffusion equation in the geodesic cylindrical coordinate, where 
longitudinal rapidity y, transverse rapidity ^ and azimuthal angle are used. The longi- 
tudinal rapidity y, and the transverse rapidity ^ are defined respectively as, 

y = \n , <^=ln 



where E,pi, pj- and m denote energy, longitudinal momentum, transverse momentum. 



and mass of the observed particle, respectively, and mj- = yVr +"^^- 
The diffusion equation in the geodesic cylindrical coordinate is given by, 

dt ~ cosh^^ dy^ + sinh^ cosh^ (^smh^ cosh^ j + ^.^^2 ^ ^^2 ' 

where D denotes a diffusion constant. However, we cannot solve Eq.© analytically at 
present. Therefore, we should consider somewhat simpler case. 

We have proposed the relativistic diffusion model, and analyzed large pj distributions 
for charged particles in Au +Au collisions 101], where a radial flow effect is not included. 

In section 2, the relativistic diffusion model is briefly explained. A distribution of 
the initial radial rapidity is taken into account. It would correspond to the distribu- 
tion of radial flow. In section 3, large pr distributions for charged particles observed 
at RHIC [3, 4] are analyzed. The magnitude of radial flow is estimated from pr distri- 
butions. In section 4, high transverse momentum limit of our model is taken and relation 
to Hagedorn's model for pr distribution inspired by the QCD is discussed. Final section 
is devoted to summary and discussions. 



DIFFUSION EQUATION WITH RADIAL SYMMETRY IN THE 
THREE DIMENSIONAL RAPIDITY SPACE 

For simplicity, we consider the diffusion equation with radial symmetry in the geodesic 
polar coordinate system, 

-T- = T- ^ smhV , (3) 

dt sinhV dp\ dp J' 

with an initial condition, 

/(p,, = 0) = ^(^. (4) 
Atz smh p 

In Eq. (|3l), p denotes the radial rapidity, which is written with energy E, momentum 
p and mass m of observed particle as. 



(5) 



Inversely, energy and momentum are written respectively as 

E = mcoshp, IpI = \J Pl'^v't ^ msinhp. 
The solution of Eq. dSj) with the initial condition Q is given |0] by 

1 



sinh-^ 



In^/AnDt sinhposinhp 
From Eq.©, the following equation is obtained; 



/(p,0 = lim /(p,po,0 = {AnDt) e 



-i/2-Dt P 



Po^o 



sinhp 



P' + Po 
ADt 



exp 



4Dti 



(6) 
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In Eq. d?!), po denotes the radial rapidity of an emission center. It would be identified 
to the radial flow rapidity. We have estimated radial flow rapidity using Eq. dTJ) |6]. 

If the distribution of po is taken into account, it would be reasonable to assume that it 
distributes randomly with dispersion Oq', 



f{poj) = {2nai)-yh-^'o/2_P^ 



■exp 



Po 



sinh Po 

Then the distribution function of radial rapidity is given by the following equation. 



Mp.t) 



f{p,poj)f{po,t)sinh Po sin OdpodOd^ 

t) 



(2;ra2^-3/2,-cT2/2_P 



e T/- — : exp 

smhp 



2a^ 



_2 



2Dt + OQ. 



(9) 
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Equation ® is obtained form Eq. ©, if 2Dt is replaced by Oj. Parameter is 
connected to the moment of po as. 



(Po') 



po/(po,0 sinh^po sinOdpoddd^ = ao +3ao. 



When p << 1, IpI = msinhp ~ mp. Then, Eq. dH) reduces to 



/(p,0^exp[-pV(2a(0 



ait) 



2Dt. 



(11) 



If we assume that Eq. (fTTl) should coincide with the Maxwell-Boltzmann distribution, 
we have an identity, kgT = mo{t)^, where kg is the Boltzmann constant. Then, Eq. (flOb 
reduces to 



a^ = kBT/m+(^-^ + ^^9 + 4{p^) 

When, Oq « 1 or Pq << 1, Eq. (fTOt is written as, 

o^^kBT/m+{p^)/3. 



(12) 



(13) 



ANALYSIS OF LARGE Pt DISTRIBUTIONS OF CHARGED 

PARTICLES 



Transverse momentum distributions of charged particles observed by STAR 10] and 
PHENIX 10] Collaborations are analyzed by Eq. Data are taken at = n/2. In 
this case, the identity, p = = , is satisfied. 

The results on the data by STAR Collaboration are shown in Fig. [l] and Table [TJ 




FIGURE 1. (a) pT distribution for Am + Am {h+ +h')/2+X aty = 10], and (b) centrality depen- 
dence of (Pq)'''^ extracted from pr distributions 



TABLE 1. Parameters on pj distributions estimated by Eq. (|9} in Am +Am 

(/z++/!-)/2+Xat3; = 0at /?]^ = 200 GeV d 



centrality 


C 




m 


y2 /n.d.f 


00-05% 


1126.3 ±51.0 


0.332± 0.003 


0.570± 0.007 


244.9/32 


05-10% 


916.2 ±47.5 


0.332± 0.004 


0.567± 0.009 


128.8/32 


10-20% 


780.0 ±36.2 


0.351 ±0.004 


0.536± 0.008 


144.2/32 


20-30% 


590.1 ±28.9 


0.360± 0.004 


0.5 19± 0.008 


98.1/32 


30-40% 


451.5 ±22.7 


0.374± 0.004 


0.492± 0.008 


95.4/32 


40-60% 


300.99 ± 14.47 


0.402± 0.004 


0.442± 0.006 


51.7/32 


60-80% 


150.91 ±7.43 


0.444± 0.004 


0.369± 0.006 


26.5/32 


PP(NSD) 


17.654 ± 1.424 


0.455± 0.006 


0.324± 0.008 


16.1/29 



Information on radial flow is extracted by the use of Eq. (fT2t under the assumption that 
the hadronization temperature is constant irrespective of centrality. The results are shown 
in Fig. lb. The mean radial flow rapidity, {pl)^/^, depends on centrality very weakly, if 
it is less than 60%. At keT = 0. 13 GeV, {pf) ~ 0.57 for centrality < 60%. 

The results on the data by FHENIX Collaboration are shown in Fig. |2land Table |2| 
Estimated value of {pi) is somewhat larger than that from STAR Collaboration at the 
same temperature. 




FIGURE 2. (a) pj distributions for Au+Au^{h++h')/2+X at y = fT], and (b) centrality 
dependence of {PqY^^ extracted from pj distributions 



TABLE 2. Parameters on pr distributions estimated by Eq. ^ in Au +Au 

{h++h-)/2 +X at y = at y/sBi = 200 GeV H 



centrality 


C 




m 


y2 /n.d.f 


00-05% 


410.00 ±42.99 


0.277± 0.007 


0.745± 0.023 


32.1/31 


00-10% 


408.50 ±41.25 


0.284± 0.007 


0.726± 0.021 


35.5/31 


10-20% 


305.21 ± 30.84 


0.293± 0.007 


0.708± 0.021 


26.4/31 


20-30% 


265.03 ±27.18 


0.313± 0.007 


0.653± 0.019 


27.8/31 


30-40% 


209.37 ±21.82 


0.330± 0.007 


0.611± 0.019 


16.8/29 


40-50% 


143.61 ± 15.42 


0.338± 0.007 


0.587± 0.018 


15.4/29 


50-60% 


103.03 ± 11.65 


0.358± 0.008 


0.537± 0.018 


9.9/29 


60-70% 


70.343 ± 8.419 


0.380± 0.009 


0.485± 0.017 


5.1/29 


70-80% 


40.688 ± 5.284 


0.390± 0.010 


0.457± 0.018 


7.4/28 


80-92% 


22.286 ± 3.126 


0.406± 0.011 


0.419± 0.018 


10.3/27 



HIGH TRANSVERSE MOMENTUM LIMIT OF THE MODEL 

In the 1970's, when large pr distributions are observed in accelerator experiments, many 
models, which have power-law behavior in pj are proposed. In Ref. [7], a model for pj 
distribution, inspired by the QCD, is proposed as {pq/{pt + Po))"- It approaches an 
exponential distribution of pr for pr 0, and a power function of pr for pr °°. 

In this section, high transverse momentum limit of Eq. Q is examined at = 
n/2, where the identity p = ln((mr -|- pT)/m) holds. Radial rapidity contained in the 
Gaussian-like part in Eq. ^ is rewritten as. 



exp 



2a| 



niT + Pt 



m 



(14) 



For any positive number e, we have, 

lim A= \ln n+P' )=0. 

pt^°° Pj pt^°° pj\m/ 

Therefore, we can approximate p in the exponent in Eq. (fT4l) as a constant, which is 
written by Icq, within some finite transverse momentum range. Then Eq. Q reduces to, 

, / m \co/o^ m ( m \co/ct2+i 
/cP,0~ ^ ■ (15) 

From Eq. (fT51) . one can see that Eq. ® shows the power-law behavior in the high 
transverse momentum limit, and that the power becomes smaller as a|, which should 
increase with the colliding energy increases. 



SUMMARY AND DISCUSSIONS 



In order to analyze large pT distributions of charged particles observed at RHIC, the 
relativistic stochastic process in the three dimensional rapidity space, which is non- 
Euclidean, is introduced. The solution is Gaussian-like in radial rapidity, where the radial 
flow rapidity po is included. It is very similar to the formula proposed in Ref. [1] at 
= :?r/2. ( See also Ref. ISfl. ) 

Transverse momentum distributions for charged particles in Au -\-Au collisions at j = 
at ^/sMN = 200 GeV are analyzed. From the observed pj distributions, an effect of radial 
flow is subtracted by the use of assumption that the hadronization temperature does not 
depend on the centrality. 

The averaged radial flow rapidity, (pQ ) depends on the centrality very weakly from 
0% up to 60%. It decreases rapidly, as the centrality increases above 60%. 

We have derived that our formula Q shows power-law behavior in pj in the high 
transverse momentum limit. Therefore, it behaves like the Gaussian distribution in pj 
when Pt « m, and like the power-law distribution in pr when pj » I- 
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